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ABSTRACT 

We calculate the reverse shock (RS) synchrotron emission in the optical and the radio 
wavelength bands from electron-positron pair enriched gamma-ray burst ejecta with 
the goal of determining the pair content of GRBs using early time observations. We 
take into account an extensive number of physical effects that influence radiation 
from the reverse-shock heated GRB ejecta. We find that optical/IR flux depends very 
weakly on the number of pairs in the ejecta, and there is no unique signature of ejecta 
pair enrichment if observations are confined to a single wavelength band. It may be 
possible to determine if the number of pairs per proton in the ejecta is >100 by using 
observations in optical and radio bands; the ratio of flux in the optical and radio at the 
peak of each respective reverse-shock light curve is dependent on the number of pairs 
per proton. We also find that over a large parameter space, RS emission is expected to 
be very weak; GRB 990123 seems to have been an exceptional burst in that only a very 
small fraction of the parameter space produces optical flashes this bright. Also, it is 
often the case that the optical flux from the forward shock is brighter than the reverse 
shock flux at deceleration. This could be another possible reason for the paucity of 
prompt optical flashes with a rapidly declining light curve at early times as was seen 
in 990123 and 021211. Some of these results are a generalization of similar results 
reported in Nakar & Piran (2004). 

Key words: gamma-rays: bursts, theory, methods: analytical - radiation mecha- 
nisms: non-thermal - shock waves 



1 INTRODUCTION 

Since t he discovery of the bright optical flash from GRB 
990123 iAkerlof et al.lll999lh many authors have considered 
this early afterglow emission as being due to the reverse 
shock passing thr ough baryonic ejecta from the explosion, as 
was predicted bv lSari fc PirarJ (I1999T) . If this early emission 
is indeed caused by the reverse shock (hereafter RS) passing 
through the ejecta, then this emission provides direct infor- 
mation about the composition and magnetic field strength 
and orientation in the material ejected from the inner en- 
gine, giving us clues about the nature of the outflow and the 
inner engine itself. Early optical afterglow data from GRBs 
990123 & 021211 have already been used to infer that the 
ejecta from these GRBs may have been fairly highly mag- 
netized baryonic ejecta JZhang et alj|2003t |Fa n et aljl2002t 
iKumar fc Panaitescull2003l ; iMcMahon et al.ll2004h . 

Only the observations of GRBs 990123 & 021211 
have exhibited this steeply falling off early optical emis- 
sion, while other bursts observed in the optical so quickly 
following the burst have exhibited either rising emission 
(030418), shallower fall-off (021004), or no emission at 



all jKehoe et al.ll200lt) . There have also been some Swift- 
detected bursts recently that have not exhibited the steeply 
fallin g off emission expected from t he RS (e.g. GRB 050319; 
iRvkoff. Schaefer. fc Quimbv 1 120051) or optical/IR emission 
varying on very short timescales, seemingly f ollowing the 
activ i ty of the centra l engine (GRB 041219a; iBlake et alJ 
2005; Vcs trand et alJ 120051) . There have also been some 
early optical upper limits from Swift of V ~ 21 magnitude 
withi n minutes after the burst (e.g. 050219a; ISchadv et al. I 
2005). This lack of observations of the optical flash begs the 
question: are GRBs 990123 and 021211 unique? What is the 
reason for the lack of optical flashes? These early upper lim- 
its on RS emission can provide information about the ejecta 
and help to answer these q uestions. 

In our previous paper (IMcMahon et alJl2004 hereafter 
MKP04), we estimated the RS emission and break frequen- 
cies at deceleration for a baryonic outflow for 10 bursts with 



1 Note however, that the sit uation in GRB 041219 a is not clear 
as the correlation reported in lVestrand et all 120051) should have 
been enhanced with a further temporal division of the optical 
signal. 
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burst parameters determined from afterglow modeling (as- 
suming deceleration time is on the order of the burst dura- 
tion, as occurs in 990123). We determined that a possible 
reason for weak RS optical emission after deceleration for 
most bursts is due to a cooling frequency below the opti- 
cal band at deceleration. Although the reverse shock emis- 
sion is expected to be bright (J?-mag ~ 10) at deceleration, 
when the cooling frequency drops below the observing band, 
no more emission from the ejecta is observed, with the ex- 
ception of off-axis e mission which falls of very quickly, as 
~ i/t( 2 +p/ 2 ) ~ i/t 3 (iKumar fc Panaitescul EoOOl . However, 
now that is is possible to quickly follow up GRBs with Swift 
UVOT, we are not seeing any bright optical flashes with this 
1/i 3 falloff. This suggests that the RS flux is suppressed for 



some reason. 



In the calculation in MKP04, we assumed the ejecta 
was purely baryonic. This may not be the case. If many e 
pairs are present in the ejecta (mixed with baryons), either 
ejected from the source or created by dissipation during the 
GRB, or if the magnetic field strength in the ejecta is very 
large (perhaps because the outflow was Poynting flux dom- 
inated), the RS emission could be suppressed significantly 
from the baryon dominated picture. Pair-enrichment of the 
ejecta causes reverse shock emission to be fainter than that 
expected from completely baryonic ejecta, and a very high 
ejecta magnetic field is likely to weaken the reverse shock as 
well; both scenarios provide alternate possibilities for sup- 
pression of reverse shock emission. It is useful then to deter- 
mine the defining characteristics of reverse shock light curves 
and spectra for a baryonic, lepton-enriched, or highly mag- 
netized ejecta, in order to distinguish between these possi- 
bilities when the Swift satellite and other robotic telescopes 
accumulate more multi-wavelength afterglow data at early 
times. 



Pair enric hment of the ejecta has been looked at by 
iLi et alJ (I2003T) and the case of highly magnetize d ejecta has 
been investigated bv lZhang fc Kobavashil J2005I) . In this pa- 
per, we take another look at pair-enrichment; we take into 
account that the reverse shock is mildly relativistic, and also 
carry out a self consistent cooling and self absorption cal- 
culation, which is necessary, since the ejecta self absorption 
frequency can be on th e orde r of or greater than the cool- 
ing frequency. ILi et alJ (I2003T) predicted a bright IR flash 
with pair-enriched ejecta, as bright as that expected in the 
optical with strictly baryonic ejecta; we find that the opti- 
cal/IR flux levels at deceleration have very little dependence 
on the number of pairs. Also, any change in spectral peak 
frequency or optical/IR flux at deceleration due to pair en- 
richment can also arise if the microphysics parameters in 
the ejecta are changed, thus creating a highly degenerate 
problem. 



We first describe our model of reverse shock emission 
and the new physics we have added to the RS emission cal- 
culation (Section 2), then discuss our results (Section 3). In 
Section 4, we summarize the differences between baryonic 
and pair-enriched ejecta. 



2 DESCRIPTION OF THE REVERSE SHOCK 
MODEL 

Here, we briefly describe our RS model, then focus on the 
improvements we have made to the calculation of emission 
from the RS. For this calculation, we use the equations for 
the ejecta dynamics and RS emission described in Kumar & 
Panaitescu 2003 (KP03), Panaitescu & Kumar 2004 (PK04), 
and Nakar & Piran 2004 (NP04) and include synchrotron 
and inverse Compton cooling. We have added several new 
components to the calculation, namely (1) the effect of syn- 
chrotron & self inverse- Compton scattering on electron cool- 
ing & synchrotron self absorption frequencies, (2) inverse 
Compton cooling of RS electrons by synchrotron flux of the 
forward shock (hereafter FS), (3) absorption of RS photons 
(particularly in the radio band) in the FS region as they 
traverse through the FS on their way to the observer, and 
(4) pair-enriched ejecta. 



2.1 The Standard Model: Dynamics & 
Synchrotron Emission in the RS 

The ejecta dynamics are determined by assuming an initial 
Lorentz factor To, the burst duration in the observer frame 
taRB, number density of the external medium no = AR~ S 
where R is the radial distance from the center of the ex- 
plosion, and isotropic equivalent energy E52 = i?/10 52 ergs. 
The Lorentz factor of the shocked ejecta with respect to the 
unshocked external medium is given by (PK04) 



r (l + 2r (n () /n ej ) 1/2 ) 



-1/2 



(1) 



where n e j = E/ (47rm p c 2 J? 2 FgA) is the comoving ejecta 
density and A is the lab frame ejecta width, taken to be 
ctdur + R/ (2ro), where tdur is the host galaxy frame burst 
duration, toHs/(l + z). When A ~ ctdur at deceleration, 
we define this case to be the non spreading ejecta case (also 
called thick ejecta), whereas when A ~ R/ (2F 2 ,), we define 
this case to be the spreading ejecta case (also called thin 
ejecta). 

When Fono/riej 3> 1/4 near shock crossing time (the 
time at which the RS reaches the rear of the ejecta shell), the 
RS is relativistic and the bulk Lorentz factor F oc _r( s_2 )/ 4 
and the ejecta radius R oc f 2 /( 4 ~ s ' (t is time in the ob- 
server frame), prior to deceleration time, for the case of non- 
spreading ejecta. In the case of spreading ejecta, Y^no/riej ~ 
0.3 for s = and ~ 0.5 for s = 2 at shock crossing time, 
so the RS is mildly relativistic. The radial distance of the 
ejecta from the center of the explosion at the time when the 
RS reaches the back of the ejecta shell (the shock crossing 
radius) is given in Table 1, as are the Lorentz factor and 
observer time at this radius. See Table 2 for the scalings of 
variables with observer time before deceleration, numerically 
determined for the mildly relativistic case. 

The reverse shock speed with respect to the unshocked 
ejecta as measured in the lab frame is (KP03) 



t n 



1.4 / rgnp 



1/2 



(2) 



where /3rs is the reverse shock speed measured in the lab 
frame and /3q is the outflow speed measured in the lab frame 
as well. The above equation is valid in all cases, from non- 
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relativistic to relativistic RS cases. It can be shown (PK04) 
that the radius at which the reverse shock has traversed the 
ejecta is within a factor of ~ 2 of the deceleration radius for 
both a wind and homogeneous external medium. 

The rate at which the ejecta electrons are swept up by 
the RS is given by 



dNe 

dR 



1.4 

ro" 



n 



1/2 



Np 
A 



(3) 



where N p = E/ (rom p c 2 ) is the total number of protons in 
the ejecta. For spreading ejecta, N e (R) oc R 3 ^ 2 for s = 
and R 6S for s = 2 (determined numerically) and for ejecta 
whose width is dominated by the GRB duration N e (R) oc 

fl (4- s )/2_ 

After deceleration (shock crossing time), the ejecta bulk 
Lorentz factor is assumed to evolve adiabatically following 
the Blandford-McKee solution, as described in PK04, as 



T e3 oc i? _(3 ~ 



)/2 (^( 2a - 3 ) z 



-1.5(4-s)/(17-4a) 



(4) 



where a is the adiabatic index, set to 4/3 for relativis- 
tic ejecta and 5/3 for non-relativistic ejecta, and A' = 
msLx(T e jCtdur, R/2F e j) is the comoving frame ejecta width. 
We take the ejecta expansion in the radial direction to be 
relativistic when the minimum electron thermal Lorentz fac- 
tor in the ejecta >10, and nonrelativistic when it is < 10. 

Synchrotron radiation from the shocked electrons is cal- 
culated assuming a power law distribution of electron energy 
with index p, i.e. dN^/d^y oc 7 _p for 7 > 7i r (where 7i r is 
the minimum electron thermal Lorentz factor averaged over 
all of the electrons as they cross the RS front), ebt, the frac- 
tion of energy in the magnetic field, and e' e , the minimum 
thermal energy given to the electrons; -fir = e' e (e p ) /m e c 2 , 
where e p — m p c 2 (T' — 1) is the minimum thermal energy 
per proton just behind the RS front, (e p ) ~ 0.5e p , and 



r' = r r (i - fop) 



(5) 



is the relative Lorentz factor of the shocked and unshocked 
ejecta. Since r'<2, we cannot approximate the evolution of 
7ir by assuming 7; r oc V'; we numerically calculate the evo- 
lution for radii near the shock crossing radius, the results 
given in Table 2. The comoving magnetic field strength in 
the ejecta is given by B' = [8ne B r n e j (e p ) (4F' + 3)] 1//2 . 

We note that ji r can be less than 1, since 7i r m e c 2 is 
defined to be the minimum thermal energy for electrons. 
We take into account that some fraction of the electrons 
are nonrelativistic and emit cyclotron radiation. We remove 
these electrons from the electron column density used for our 
synchrotron calculation, i.e. n co ;, r = A r e (7i r ) p_1 /47rJ? 2 , and 
set the minimum thermal energy for electrons to be m e c 2 . 

After deceleration, for adiabatically cooling electrons, 



(R 2 A') 



-1/3 



where V' is the comoving 



7*' 

ejecta volume, and the ratio of specific heats is a e = 4/3. 
The magnetic field is assumed to be predominantly trans- 
verse and frozen into the ejecta, decaying as B' oc (RA') 1 . 
When radiative cooling becomes less important than adia- 
batic cooling, the Lorentz factor of electrons cooling on a 
dynamical timescale, 7 cr (discussed in more detail below), 
evolves in the same manner as "fi r - 

The observer frame synchrotron injection frequency is 
calculated with 



0.98gB'7 2 r r 
27rm e c(l + z) ' 



(6) 



where we use here and elsewhere the common notation of 
v x being the synchrotron frequency of an electron with a 
Lorentz factor 7^. The synchrotron flux at the peak of the 
/„ spectrum is given by 



AUVXi + z) 

4nd 2 L 



(7) 



where N e = N e (R) is the number of electrons heated by the 
RS determined from Equation^] P„ — 1.04g 3 B/m e c 2 is the 
comoving power radiated per electron per unit frequency at 
the peak of /„, and d L = 2c % /TTz[(l + z) 1/2 - l]/Ho is the 
luminosity distance. We use Ho = 65 km s _1 Mpc -1 and 
for simplicity fl — 1, A = 0. The f actors of 0.98 in Equation 
HJand 1.04 above for P Vp are from lWiiers fc Galamal dl999l) 
for the case of p = 2. 



2.2 Inverse Compton Cooling & Synchrotron Self 
Absorption 

We include the inverse Compton process to calculate elec- 
tron cooling. The Compton Y parameter is obtained by solv- 
ing the equation describing radiative loss of energy for a 
single electron with a Lorentz factor y e 



d (7 e m e c 2 ) 



otcB' 2 ^ 



dt' 



67T 

fv'>u' a 
J total 



2a2 
r es r y p P p 



(1+Tsa (!/')) 



(8) 



where primes denote variables measured in the rest frame of 
the shocked fluid. The bracketed terms on the RHS are effec- 
tively (1 + Y); the first part, the fraction of energy emitted 
from one electron with Lorentz factor 7 C that is not absorbed 
by material in emitting region, is given by 

f, , 

Jv'>u'„ 



f 

J total 



£,x 4 / 3 -25S3 

3[2e" a: ^+^(l-Erf(yS))] 

e- a =(2o:3/ 2 +6v^)+3v^ : (l-Erf( v ^)) 



for v' c > v' a 
for v' c < v' a 



(9) 



where x = v' a lv' c and Erf is the error function. This equation 
is necessary for calculating electron cooling, since v' c can be 
less than the self absorption frequency, v' a , in the RS. 

The second part of the bracketed term in Equation |H1 
is equivalent to Y; r es is the Thomson optical depth to 
electron scattering, 7 P = min(7i,7 c ) and (3 P is the corre- 
sponding velocity in units of c, v' p — min(i^, u' c ), okn is an 
approximation of the Klein-Nishina correction to the elec- 
tron scattering cross section, r sa (y 1 ) is the optical depth to 
synchrotron self absorption, and / (v 1 ) is the normalized set 
of broken power-laws describing the synchrotron spectrum 
(i.e./ (v' p ) = 1). 

Both of the bracketed terms of the effective (1 + Y) 
(Equation |HJ are dependent on u' c and u' a , and in turn, v' c 
is also dependent on (1 + Y), so the equations for the three 
variables must be solved simultaneously; this has been done 
for all numerical results in this paper. It is, however, possible 
to semi-analytically estimate v' c when we assume that Y <S 1 



Table 1. Values of Parameters at Deceleration 



= 



: Ctdu 
= 2 



A = R+/2TI 
s = 



A = R+/2TI 
s = 2 



R+ 

r+ 

*obs,4- 
( e p-. + ) . 



£ l/2 
Br, 
1/2 

Br-,-1 



1 + 
-1/2 



£^(^1+1 
1/2 /2 



ir,4- ' 



X 



C 2 (l+*) 

■cr,+ =e 3 J*_ 1 (l + z)C% 2 X 

_ c /1.5 c 0.8: 
~ e er,-l e Br, 



/-1.38 /--0.5 v 



55 (n0m p 5, tL J 



1/2 



2M dur 
T (0.18nom p c 5 7rt5 ur i 
(o.91n m p c 5 7rt3 ur roE- 



.A V8 
.A 1/8 



21.5 (n 3 E 52 t-l) 1/8 G 



51 [ ^ | 4 

T (0.79A^m p c 3 t d . ar lrE- 1 ) 1, ' ^ 
(0.99Am p c 3 TTt dar r o E- i y /4 ' 
2.2 xlO* (AfB^t-^^G 



5.1 x 10 n nJ /2 r2 2 Hz 

i-exio"^ 1 ^ 2 ^ 2 ! 

12.9n«-«3 B i 2 3 r o ; 5 td - u 3/4 



5.4 x lO^r^t-^B-^Hz 
3.6 x 109^ 2 /2 t^ 2 r A- 2 Hz 

1.2 x loMl-srgJ^t^ 



0.69r 

y ri c j m p 7rr^ y 
0.03 
1.07 



1/3 



1/3 



b 



x 791.0nJ /4 ^ 2 /4 t du 3 r /4 r^ 2 mJy 1.3 x lO^^^^r^mJy 26.1n /2 ,E52ro ! 2mJy 



[7.9^™^ TrJJ 
0.83r 



2.4 x 10 10 n,!, /2 r2 j2 Hz 

2.2 x io 16 r 4 / 2 3 A--V6 B52 2 / 3 Hz 

0.043njj' 88 r2j£; 52 



0.Q2E 



0.01 
1.02 

3.3 x io 5 A 3/2 r 3 > 2 £;- 1 G 

2.6 x 10 6 A 3 / 2 r 3 1 2 mJy 

2.6 x lO^A'^rj^B^Hz 
5.3 x 10 8 S 5 2A^ 5/2 r- 2 Hz 
4.6 x W 6 A 2 , 6 rl 2 E-°- 75 



A, = A/10 35 , E 52 = -E/10 52 , r , 2 = To/100, ! = e' er /0.1 and e Br ,-l = e Br /0.1. f VyR is the flux at deceleration in mjy in the R-band for p = 2.5 and z 

{fir, v a r) < vr < v cr is assumed. For full expression used in Figure 5, see Appendix A. 
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Table 2. Scalings Before Deceleration 





A ctdur 


^ ddur 


A = R+/2Y 2 


A = R+/2Y 2 




s = 


s = 2 


s = 


s = 2 


R 


t l/2 


t 


t 0.67 


t 0.88 


r 


t -l/4 


t° 


t -0.17 


j-0.066 


7ir 


t 0.4B 


t° 




t 0.62 


fv,X 


t -0.18+0.19p 


t l-0.5p 


t 1.5(p-l) 


t 0. 012+0. 082p 


fv,R 


t 0.29+0.19p 


t -0.5(p-l) 


£-0.86+1. 5p 


j-0.56+0.082p 



Values of exponents in this table are determined numerically for 
mildly relativistic spreading ejecta and determined analytically 
for relativistic non spreading ejecta, near, but prior to, deceler- 
ation time, f is time measured in the observer frame. For the 
non-relativistic case of Y 2 n /n eJ « 1/4, T a Ji° and R oc t. 
fv.x is the synchrotron flux for v bs > (fin^cr^nr), and f v R 
is the synchrotron flux for (ui r ,u ar ) < v B j, a < v cr . Scalings after 
deceleration can be found in e.g. PK04. Note: it is possible for f u x 
to be decaying before deceleration because it peaks well before 
deceleration. 



and IC does not contribute to the cooling calculation (see 
Table 1). 

The comoving frame synchrotron self absorption fre- 
quency is calculated by equating the comoving synchrotron 
flux at v' a to flux from a black body in the Rayleigh- Jeans 
part of the spectrum, or 



2kT- 



(10) 



where kT — max [y a , min(7;, 7 C )] m e c 2 . A frequent arrange- 
ment of the RS break frequencies at deceleration is v[ < 
v' a < v' c ; in this case, 



JL 

2m a 



( qB' 

\ 2itm e 



1/2 



(11) 



where f p = N e VSq 3 B'/ (inR 2 m e c 2 ^j is the comoving flux at 
the peak of f v . 

The scalings of RS optical and X-ray flux with observer 
time just before deceleration for the four cases of s = 0, 2 and 
A = (ctdur, R/2Tq) are shown in Table 2. Only synchrotron 
emission is included for these scalings (inverse Compton is 
also important in the X-ray), and Y <C 1 is assumed (the 
scaling for f Ut R is not affected, but f U} x may be). 

After deceleration, if ejecta is in the radiative regime, 
we continue to calculate the RS cooling and synchrotron self 
absorption frequencies v cr and v ar and the effective (1 + Y) 
by solving Equations Icl fc HTH simultaneously, as done before 
deceleration. If adiabatic cooling is more efficient, y cr decays 
as ji r (§ 12. H and the self absorption frequency is solved for 
using Equation llUI After v cr falls below the observing band, 
the ejecta emission turns off, and the observed flux is due 
to off axis emission that de cays approximately as f~( 4 + p )/ 2 
jKumar fc Panaitescdl200oTl . 



2.3 Inverse Compton cooling by an external 
source 

We have also allowed for external sources of flux to influ- 
ence the cooling of electrons in the ejecta. In particular, we 



include synchrotron flux from the FS incident on the ejecta 
in our calculation of v cr . The external synchrotron flux in 
the shocked ejecta comoving frame is given by 



/ex — 



t rft' 2 <v 2 



(:)7T 



(12) 



where the subscript "ex" denotes values from a source ex- 
ternal to the cooling calculation being done (here, FS flux 
is the external flux incident on the RS ejecta), "f 2 i£X is the 
average electron thermal electron Lorentz factor squared in 
the forward shock. The equation for 7 C which includes the 
contribution of f ex to the electron cooling is: 



Xi (1 + 50 + 



where xi — arB' 2 t' / (67rm e c 2 ), ^2 = crrt' / (m e c 2 ) , and 
o~KN,ex = [hvp tex "/ c / (m e c 2 )] 1+Qex i s an approximation to 
the Klein-Nishina correction on the external flux, and u P:ex 
is the peak of the vf v spectrum in the FS; a ex = —0.5 if 
Vi,ex > fc,ei, and (1 — p)/2 if i/ i>ex > u C)EX . We also allow for 
FS synchrotron self-Compton flux to influence RS cooling; 
however the high energy (3> m e c 2 ) of the FS synchrotron 
flux suppresses self-inverse Compton scattering, thus not sig- 
nificantly contributing to the RS cooling calculation. We do 
not consider the effect of RS emission on FS cooling; at most 
it is an order unity effect. 

2.4 Absorption of RS Photons in the FS Material 

We include the effect of absorption of RS photons in the 
FS material. For a selected observing band, v ^ B , the optical 
depth to absorption is 



T a bs,FS = < 



5/3 



v b a < min(f /, Vif) 

f \-(p+4)/2 (I 4 ) 

\-^fj v °bs > min(f c /,^/) 

where v a f is the observer frame FS synchrotron self absorp- 
tion frequency. Before deceleration, the shocked ejecta and 
FS medium are moving together; however, after decelera- 
tion, these regions are moving at different Lorentz factors, 
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and one must take care to use the appropriate value of p a f 
when the RS photons are passing through the FS. 

After calculating the RS emission as described above, 
the absorption in the FS is taken into account by reducing 
the RS flux by a factor of exp(— T a b a ,Fs)- This exponential 
cut off in flux can significantly reduce the RS emission, es- 
pecially for observations at longer wavelengths, such as the 
radio. 

2.5 Lepton-enriched ejecta 

We include e pairs in the calculation by simply adding a 
certain number of pairs per ejecta proton, N±, to the calcu- 
lation described above. We do not calculate the creation of 
these pairs; we make the assumption that a certain number 
of pairs per proton are present in the ejecta already, either 
being intrinsic from the source or being created by dissipa- 
tive processes during the GRB, prior to the afterglow stage. 

We change the radiation calculation described above by 
multiplying the column density of electrons by the number 
of ejecta pairs, L = 1 + 2N± (the 1 is for the electrons 
already accompanying the protons in the baryonic ejecta) 
and by dividing the minimum electron energy by the num- 
ber of pairs (i.e. the new minimum thermal Lorentz factor 
is 7i r /L). The dynamics calculation is only altered when 
N±>m p /m e . To account for the presence of a high number 
of pairs, we alter the dynamics calculation by reducing the 
number of protons in the ejecta for a fixed burst energy by 
L ep = 1 + 2N±m e /m p (see Table 1). 

The RS minimum electron thermal Lorentz factor av- 
eraged over the entire lepton population at shock crossing 
is 



Since the RS is only mildly relativistic in the spreading 
ejecta width case, the addition of pairs to the ejecta quickly 
drops 7i r into the Newtonian regime. The already low in- 
jection frequency in the reverse shock can be reduced dra- 
matically by even modest L (see Table 1). The injection fre- 
quency of course cannot drop below the cyclotron frequency; 
we keep track of this in our numerical calculation. 

The cooling frequency has very weak dependence on 
N± (there is dependence on N± through L ep , but this only 
makes a difference if ./V±>1000). If y > 1 and Vi r < v cr , 
v cr oc i*(vr-2)/{A- Pr ) ( for N± < mp /m e ). The depen- 
dence of p ci r on the pair content when Vi T < p ar < Per is 
Par oc ~L,( 4 ~ 2 Pr)/(4,+Pr) f or jy ± m p /m e ; Par is also fairly 
insensitive to pair content. 



3 RESULTS 

In this section, we describe the effects of IC cooling by flux 
generated in the RS and FS, absorption in the FS, and the 
effect of lepton-enriched ejecta on the RS emission, then de- 
termine if there are observable signatures of pair enriched 
ejecta. To ascertain the effect of each of these new addi- 
tions to our calculation over the entire parameter space and 
a wide range of RS strengths for s — 0, 2, we randomly 
vary each parameter in the ranges: 50 ^ To ^ 1000, Is ^ 
t G RB < 100s, 0.1 < E 52 s? 1000, 10" 3 cm -3 < n ^ 100 




0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 
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Figure 1. Left panel: Ratio at deceleration of numerically calcu- 
lated Per (with IC cooling) to analytic estimate given in Table 1 
(not including IC cooling) for 1000 test cases in s = 0. The inclu- 
sion of IC cooling reduces p cr by a factor of roughly 10 over a wide 
range of RS strengths (£). Right panel: Ratio at deceleration of 
numerically calculated cooling frequency without external cool- 
ing included to numerical value with external cooling included, 
for the also for the s = case. The legend gives the value of 
Egf in comparison to the value of e^ r , e.g. for the red squares, 
?Bf = O.leBr- 

cm" 3 (10~ 2 < A, < 10 for s = 2), 10~ 5 sC s Br < 0.1, 
0.01 < e ir < 0.1, 2.01 sC p r < 2.91, and 1 sC N± «S 1000 for 
1000 test cases (synthetic GRB afterglows). All of the pa- 
rameters are assumed to have uniform distributions in log 
space with the exception of p r , which is assumed to have 
a uniform distribution. The microphysics parameter ranges 
were chosen to be consistent with values found from late 
time afterglow modeling. 

3.1 Inverse Compton & External Cooling 

In the left panel of Figure^ the ratio of our numerical value 
of p cr ,+ to the analytical value given in Table 1 for 1000 test 
cases is plotted against £, the dimensionless RS strength pa- 
rameter, at deceleration, for s = 0. 2 Each RS parameter is 
randomly varied in the ranges described above, except for 
iV±, which is set to 0. For each of the 1000 cases, the spread- 
ing condition is evaluated at deceleration, and the numeri- 
cally calculated p cr ,+ is compared to the analytical value in 
the proper ejecta evolution case. 

The numerically calculated value of p C r,+ , which in- 
cludes the self consistent inverse Compton and synchrotron 
self absorption calculation as described in H2.2I is in general 
about a factor of 10 lower than the analytic value (not in- 
cluding IC losses) that is typically used over a wide range 
of RS strengths. For those cases where the numerical value 

2 £ = (l/A) 1 / 2 r~ 4,/3 , where I = (3£/(47rn m p c 2 )) 1 / 3 ; £ < 1 is 
relativistic RS and £<1 is the generic, mildly relativistic RS case. 
For £ ru 2, the ejecta is in the spreading regime with a mildly to 
non-relativistic RS (these points are bunched very closely together 
in the figures, since £ is virtually constant over the dynamics pa- 
rameter space for all spreading cases), whereas the non spreading 
ejecta case ranges from { < 1 to £5;2. 
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is even smaller, the Compton Y parameter is rather large 
(these cases are more abundant for the non spreading case, 
with relativistic RS). For s = 2 (not shown in the figure), the 
analytical value of v cr ,+ can also be less than the numeri- 
cal value of v cr ,+ . This is because our cooling calculation 
has taken into account the fraction of the synchrotron flux 
that has been synchrotron self-absorbed, reducing the rate 
of energy loss for electrons; synchrotron self-absorption con- 
tributes more to the RS cooling calculation in s — 2 than 
s = 0. 

External flux influences the RS cooling calculation by 
decreasing v ar . If the FS radiation was produced with sbj ~ 
£Br, the external flux can lower v cr via IC scattering by up 
to a few orders of magnitude for both s — 0, 2 cases. The 
effect of external flux IC cooling on the RS synchrotron flux 
near deceleration in most cases is relatively small; the effect 
is greatest in the x-ray band, where v bs > (y a r,Vi r ,v cr ), 
and flux here can be decreased by a factor of a few (IC 
flux is also an important contributing factor to the total 
light curve in the x-ray, and is not considered here). The 
more important effect of the external flux is to decrease the 
cooling frequency below the R band on a shorter time scale 
after deceleration, causing the light curve to fall off very 
steeply, as ~ 1/t 3 , soon after deceleration. Overall, the effect 
of external flux on electron cooling makes RS synchrotron 
emission a little fainter at the peak of the light curve and 
fall off faster after deceleration, making the RS synchrotron 
emission more difficult to detect. 

In the right panel of Figure Q we show the effect of 
adding external cooling to our calculation. We plot the ratio 
of the numerically calculated v cr ,+ without external cooling 
to v c r,+ including external cooling from the FS. We let pf = 
p r , Sif — Sir, and try three cases of e B f = (1, 0.1, 0.01)esr- 
With external cooling included, v cr in either the spreading 
or non spreading case can be reduced by up to 3-4 orders of 
magnitude! This is true over a wide range of RS strengths, 
and also for s = 2. The effect is strongest when es/ = EBr 
and weakest when £s/ = 0.0l£B r . 



3.2 Effect of Absorption in FS 

Absorption of RS synchrotron photons in the FS is most 
important in the radio, as the FS self absorption frequency 
lies in this band. In a large part of the parameter space, 
the low frequency (~4 GHz) RS flux can be completely ab- 
sorbed; in a another run of 1000 different RS cases (varying 
the parameters as described above, and setting FS Ej/, £ B f, 
and pf equal to RS values), we find that in approximately 
15% of 1000 test cases, the RS radio (8.5 GHz) flux is com- 
pletely absorbed by the FS material. In approximately 40% 
of the test cases, r abS: FS > 1 (including those cases where 
the RS radio flux is completely absorbed) , meaning that the 
RS radio flux was decreased by at least a factor of 3. If we 
set EBf — esr/100, we find that 12% of 1000 test cases are 
completely absorbed in the FS, and in 23% of the cases, 
T a bs,FS > 1- Absorption in the FS may turn out to be a 
contributing factor to the difficulty in observing RS flux in 
the radio. 



3.3 RS Emission with Pair-Enriched Ejecta 

Adding a certain number of pairs per proton to the ejecta 
adds another parameter to the RS flux calculation, bringing 
the total number of parameters in the RS to four, viz. Ei r , 
EBr, Pr, and N± 3 . We need four independent measurements 
of the RS emission to constrain these four parameters. Flux 
in the optical & X-ray at early times, spectral index and ra- 
dio flux at ~ 0.1 day are four quantities that can be observed 
with the current generation of instruments. 

For v obs > {vir,v cr ,v a r) (x-ray) with Y < 1 or v ir < 
Var < Vobs < Vcr (optical) for any Y, the observed flux at 
deceleration has a dependence on N± of 

L (2-Pr) 

f" K T 5(2- Pr) /4 - ( 16 ) 

For p r ~ 2, this flux depends very weakly on the number 
of pairs in the ejecta. For p r > 2, the x-ray and optical 
flux are decreased from the N± = case as N± is increased 
(for N± <C m p /m e ), but the shape of the light curve does 
not change (i.e. no breaks are introduced by reducing the 
injection frequency below the low value at which it already 
sits in the absence pairs in the ejecta). 

In the radio, the observer flux for v obs < u ir < v ar < v cr 
is proportional to L _1 L,;p. The flux in the radio is the most 
greatly affected by the addition of pairs to the calculation, 
so the best possible place to look for a signature of pair- 
enriched ejecta is in the radio wavelengths. In Figure [5] wc 
plot an example of RS and FS radio light curves (8.5 GHz) 
without pairs (for burst parameters see figure caption) com- 
pared to three cases: (1) adding 100 pairs/p to the ejecta, 
(2) decreasing EBr by a factor of 10, and (3) decreasing Ei r 
by a factor of 10. All three of these actions produce very 
similar RS radio light curves; the situation is highly degen- 
erate, and it would be difficult to distinguish the effect of 
pairs in the ejecta on the emission from a variation in the 
microphysics parameters in the RS. 

1000 RS test cases (parameters varied over ranges 
described above with no external cooling or absorption 
in the FS) show that the peak frequency of f v , at 
ma,x{v ar , min(vi r , v cr )}, at deceleration also has very little, 
if any, dependence on the number of pairs in the ejecta for 
s — 0, 2. This is also true if external cooling is included in 
the calculation. Th is is in contradiction to the conclusion 
that iLi et alJ i200Sf) came to; this is because, over a large 
portion of the parameter space, the peak of the spectrum in 
both the baryonic and pair-enriched cases is at v ar , and v ar 
has little dependence on the number of pairs in the ejecta. 

Thus, we conclude there is no single wavelength light 
curve signature that depends only on the pair content of the 
ejecta, and not on the microphysics parameters in the RS as 
well. 



FS microphysics parameters and the energy in the explosion, 
external density, and initial jet opening angle can be determined, 
as has been done in the past, by late-time FS light curve fitting 
when the RS makes a small contribution to the total flux. 
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Figure 2. Radio (8.5 GHz) light curve for s = with the input 
parameters £52 = 3, To = 160, tQRB = 4 s, no = 1.6 X 10~ a 
cm -3 , e Br = e B f = 4 x 10 -5 , £j r = Sif = 0.09, p r = Pf = 2.5, 
N± = 0, and 2 = 1. External cooling and absorption in the FS 
are not included in this calculation. The solid line represents the 
RS light curve with the above parameters and the FS emission 
is shown by the dotted line. The three colored lines alter the 
parameters in the calculation as detailed in the legend. 



3.4 Using Optical and Radio Observations to 
Determine Pair Content 

The best way to discriminate between baryonic and pair- 
enriched ejecta is by comparing flux at low frequencies (radio 
band) and the optical band. NP04 suggested a particular 
combination of the ratio of observed RS optical and radio 
fluxes at the peak of each respective light curve and the ratio 
of the observer times at which these peaks occur as a way of 
determining if the radiation is from the reverse shock. We 
use a similar combination to decide whether the ejecta is 
pair enriched. 

The parameter we use to determine the pair-enrichment 
(x) is slightly different than that in NP04 (note the differ- 
ence in exponents on the ratio of observer times), 

* S (t) = (t^)^ ~ 1000 (17) 

\r0/ \toJ \VradioJ 

where 



5( Pr 



ill + 2 



(Pr 



s = 
s = 2 



(18) 



and are determined analytically using decay indices after 
deceleration of Vi r oc t~ J ^ z and F pr oc t~ 9 ^ s for s = and 
Vi r oc t~ 2 and F pr oc t~ 5 ^ 4 for s = 2 (assuming spread- 
ing ejecta evolution for both cases). As in NP04, to is the 
observer time peak of the optical RS emission, t* is the ob- 
server time peak of the radio RS emission (the peak is pro- 
duced when v ar falls below the radio), and Fo & F* are the 
observed fluxes at these two points. The shell width at de- 
celeration in both shell evolution cases is ~ R+/T+ which is 
<C R 4 , so time delay between photons arising from the front 

4 This is true in the non spreading case because of causality-the 
distance traveled by the shock front in the shell comoving frame 
is of order i?+/r+, which is iJ+/T? in the lab frame 



Figure 3. \ f° r Pr = 2.5 and s = for 1000 test RS numerical 
light curves, excluding absorption of RS photons in the FS and 
external cooling. We randomly vary all RS parameters but p r 
and N± ; other values of p r give qualitatively similar results when 
adding leptons to baryonic ejecta, however \ for the reference 
purely baryonic case is higher for higher p r and vice versa, s = 2 
has qualitatively similar results which are discussed in i|3.4l The 
scatter in each case is caused for the most part by the variation in 
dynamics parameters. There are some points with very low \ cut 
off of the bottom of the plot; these are caused when v ar > v cr , 
and v ar tracks Vcr in time (which evolves more quickly, leading 
to smaller i„). 



and back ends of the shell has small effect on the observed 
flux. 

We plot x i n Figure [3] versus the dimensionless RS 
strength parameter, £, in a run of 1000 RS cases as described 
in § |21 except that p r is held fixed at 2.5 and N± is held 
constant for each case. External cooling and absorption in 
the FS were turned off. For N± = 0, the value for \ can De 
larger or smaller by a factor of a few from 1000 for a constant 
value of p r ~ 2.5, and is fairly constant over a large range 
of reverse shock strengths. There are a few points which are 
cut off the bottom of the plot with abnormally low x~th ese 
points are caused by v cr falling below v ar . When this hap- 
pens, v ar tracks Vc r (since there are no electrons radiating 
above v cr after the RS reaches the rear of the shell), which 
evolves more quickly than v ar does when v ar < v cr , making 
i» and F* smaller. This scenario is easily detected by check- 
ing the decay of the radio light curve after the peak; if the 
decay is steep, like ~ 1/t 3 , the abnormally low value of x is 
due to v cr falling below the radio band with v ar - 

When pairs are added to the ejecta, we find that for 
iV±>100, the value of x has dropped by a factor of about 
5 for the parameter space with £<1.6 (ultra-relativistic to 
mildly relativistic RS). For N± = 10, the value for x ranges 
from ~ 3000 when the RS is relativistic (£ < 1) to ~ 200 
when the RS is mildly to non relativistic, making it difficult 
to distinguish between baryonic and pair enriched ejecta. 
Compared with the range of y for baryonic ejecta (greater 
or less than 2000 by a factor of about 2), the drop in y by 
a factor of 5 for iV±>100 could possibly be used as a tool to 
distinguish the pair content of the ejecta. So, if observations 
are available for a burst at the peak of the R band and 
radio RS emission, and p r is able to be determined from 
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Figure 4. Left panel: Ratio of \ (no pairs) without external flux 
or absorption in the FS to \ with external flux and absorption in 
the FS plotted against the value of £ at deceleration. Note that 
the ratio including absorption can go to arbitrarily high numbers, 
since radio flux can be completely absorbed (x = with absorp- 
tion in FS in these cases); some of these points have been cut off 
the plot. Right panel: Ratio of \ calculated including contribu- 
tions of FS in just the radio and in both the radio and optical 
band to \ calculated without FS contribution, plotted against £ 
at deceleration. 



spectra, we can calculate \ an d determine if iV±>100 (it 
is not possible to determine the exact number of pairs per 
proton, but only if iV±>100). It will be extremely difficult 
to tell if the ejecta has a pair content of iV±<100. 

In Figure 2] in the left panel, x is shown including the 
effect of adding external cooling and FS absorption to the 
purely baryonic ejecta case. We repeat the calculation of 
1000 test cases with p r = 2.5 and N± = 0, and we set 
Pf — p r , and vary and ebj in the same ranges as the cor- 
responding RS parameters. The inclusion of external cooling 
in the calculation does not affect \ much. The value of x is 
increased at most by a factor of a few from the calculation 
without external cooling included, and still lies within the 
scatter for the baryonic case shown in Figure [J] 

Including absorption in the FS, however, will invalidate 
Equation 1171 since we assumed in deriving it that the RS 
radio flux at the light curve peak is equal to the flux at the 
peak of the RS spectrum at t, . This decrease in F* by a large 
factor also decreases x by the the same factor (x in this case 
could be <1). In >I3.2I we found that absorption in the FS 
can be a large effect for much of the parameter space at 8.5 
GHz. Absorption in the FS for observer frequencies less than 
~ 10 GHz may indeed make it more difficult to accurately 
determine the pair content of the ejecta. However, at higher 
frequencies, ~ 100 GHz, FS absorption is less important 
and one can calculate x a t these frequencies to determine 
N± (note that x will be a factor of a few smaller at higher 
radio frequencies for the purely baryonic case). 

To determine if x for s — 2 has similar properties to the 
s = case, we carry out a numerical simulation of 1000 RS 
test cases for s = 2 identical to that done for s = (without 
external cooling or absorption) , finding qualitatively similar 
results as the s = case. With iV±<100, the pair content 
is difficult to determine, but for iV±>100, the value of x is 



about a factor of 5 lower than in the baryonic case. In the 
purely baryonic case, x ~ 3000 for s = 2 and this value 
ranges between a factor of 3 higher and lower. The spread 
of x values here is larger because the range of no = AR~ S 
for s — 2 is larger than the range of no chosen for the s = 
case. 

Up to this point, we have not considered the contami- 
nating effect of the FS emission to the optical or radio RS 
emission at the RS light curve peaks. It may be most diffi- 
cult to distinguish RS from FS flux in observed radio light 
curves, so we calculate the value of \ with the FS radio flux 
included. In the right panel of Figure 0] we've plotted the 
ratio of x calculated (for s — 0) with the FS radio flux con- 
tribution to x without FS flux at the time which the RS 
radio light curve peaks. Again, for 1000 test cases, we set 
Pf = p r , set N± — 0, and vary Eif and £s/ in the same 
ranges as the corresponding RS parameters. We find that 
the value of x can increase by a factor of up to 100 from the 
RS only case, outside of the scatter in x- I n 60% of the 1000 
test cases, the FS radio contribution was greater than the 
RS radio contribution at the RS radio peak. For s = 2, this 
occurs in 83% of the 1000 test cases. One needs to separate 
the contribution of the RS and FS to the radio flux using 
late time data in order to use this tool reliably to determine 
the ejecta pair content. 

If the RS radio peak was observed, but the light curve 
has significant FS contribution, it may be possible to sepa- 
rate the contribution of FS from the RS radio flux by con- 
tinued monitoring of the radio band for a period of a week 
or so, when the FS radio LC peaks. This information can be 
used to determine the contribution of the FS to the radio 
flux at the time of the RS peak, and x due to the RS alone. 
If the radio flux is dominated by the FS contribution during 
the RS peak and the RS peak is not observed, then this tool 
cannot be used to determine ejecta pair content. 

We also find that the optical RS peak may be as diffi- 
cult to observe as the radio peak due to FS contamination, 
in agreement with the conclusion drawn in NP04. In the 
right panel of Figure 0] we also show the ratio of x with FS 
contribution (both radio and optical) to x with only RS flux. 
The inclusion of FS optical flux at deceleration increases the 
scatter in x even more; % with FS flux ranges from being 
100 times smaller to 100 times larger larger than x with 
RS flux only. For the s — case, the FS optical flux at de- 
celeration is brighter than the RS optical flux 70% of 1000 
test cases; in the s = 2 case, this occurred in only 30% of 
1000 test cases. This is another difficulty in calculating x 
for observed light curves. However, as in the radio, if an RS 
optical peak is observed, the FS contribution may be able 
to be removed if the optical light curve is followed for long 
enough after deceleration to determine the FS contribution 
to the total observed flux. We note that the fact that the 
FS optical emission is larger than the RS optical emission 
at deceleration in ~70% of test cases for s = may explain 
the lack of optical flashes and rapidly declining light curves 
at detected early times. 

Scintillation in the radio may also be a problem for ob- 
serving the RS radio peak due to the high level of vari- 
ability that this process introduces into the obs erved radio 
light curve, as was obse rved in GRB 970508 dFrail et alJ 
Il997t iTavlor et al.lll997l) . The fluctuations are more pro- 
nounced at early times, when the RS radio light curve is 
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expected to peak, because the source size is smaller. One 
can reduce fluctuations from scintillation by observing at 
frequencies higher than 8.5 GHz. Observations made at fre- 
quencies near 50 GHz and even in the millimeter range (250 
GHz) may be more suitable for determining \ an d N± be- 
cause of the insensitivity to scintillation at these frequencies. 
The effect of absorption in the FS is also smaller for these 
higher energy radio photons (T a bs,FS > 1 for a much smaller 
parameter space), and the RS light curve peaks at earlier 
times at these frequencies when the FS light curve may not 
be as bright. With current technology, radio observations 
on the timescales necessary for calculating \ °f f ew min- 
utes/hours to days after the burst are feasible at around 
8.5 GHz. At higher frequencies like the sub-millimeter, re- 
cent observ ations are typically being done ~ 0.1 to 1 day 
post-burst JSmith et a,lJl2005T : the timescale for the sub- 
millimeter light curve peak is from a few minutes to a few 
hours after the burst (up to about 0.1 day). 

Other things that may limit the usefulness of \ as a 
tool for estimating N± include sparseness of sampling of the 
radio light curve near the RS and FS peaks and breaks in the 
electron energy distribution between the optical and radio. 
The sparseness of radio data points near the RS and FS light 
curve peaks may introdu ce an error into the determin ation 
of x °f a factor of a few. IPanaitescu fc Kumarl J2002) have 
found evidence in late-time afterglow modeling of a break in 
the electron energy distribution between the radio and the 
optical; one must be careful to choose the correct value of 
p r (value of p r between vt r and the optical band) for the 
calculation of X- 

In summary, it is possible to use optical and radio flux 
and observer times at the respective peaks of the light curves 
to determine pair content if the number of pairs/proton is 
>100. There are many factors, however, that one must take 
into account when using this tool. The contributions to the 
optical and radio flux from the RS and FS must be sepa- 
rated, which may be difficult and requires good time cover- 
age from ~ 1 minute to hours in the optical and ~ 1 hour to 
days in the radio. In the radio band, below 10 GHz, absorp- 
tion in the FS and interstellar scintillation pose problems, 
and in order to avoid these issues, one should use observa- 
tions at a higher frequency than ~ 100 GHz. 

3.5 Constraining Parameters with Available 
Observations and Upper Limits 

In Figure we plot the comparison of our numerically cal- 
culated R band flux at deceleration with the full analytical 
expression given in the last row of Table 1 (p = 2.5 and 
z — 1) and in Appendix A (full expression). The numerical 
calculation does not include external cooling or absorption 
in the FS (the latter is unimportant for the R band, but 
the former could decrease the flux by a factor of ~ 2), and 
all parameters are varied, including N±. The analytical ex- 
pressions assume (v ir ,v ar ) < v b s < v CT and include pairs. 
For each case of s — or s — 2, we ran 1000 test cases, 
and separated these 1000 cases into two sets determined by 
their ejecta width at deceleration. The blue line shown on 
each plot shows a linear relationship for a guide, and the 
orange line shows the relationship if the numerical flux were 
10 times the analytically calculated flux. 

In all four cases, the scatter around the linear relation- 
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Figure 5. Comparison of 1000 test cases of numerically calcu- 
lated R-band flux at deceleration with analytical expression for 
each case of s = 0, 2 spreading/non-spreading. All parameters, 
including N±, are varied. The full analytical expressions assume 
(ui r , v a r) < I'obs < v cr and are given in the last row of Table 1 for 
the case of p = 2.5, 2 = 1 and in Appendix A for the general case. 
A linear relationship (blue line) and fn analytical =10 times fn 
numerical (orange line) are shown on each panel for a guide. 
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Figure 6. Observer frame RS flux in the ij-band at de- 
celeration assuming z = 1 and including pairs (no external 
cooling or absorption in the FS) plotted against the param- 
eter dependences in the expression for flux in the R band 
at deceleration (without constants). All parameters are var- 
ied. The variable D in the expression on the x-axis contains 
all of the dynamics parameters and is different for each of 
the four cases shown: D = n, 
iPr/2 
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dashed lines denote the observed R band flux at deceleration for 
G RB 990123 shifted to 2 = 1, observed upper limits reported 
in iKehoe et a7] <200ll) from the ROTSE telescope, and the V- 
band upper limit re ported by Swift UVOT for GRB 050219a 
( Schad v et al. shown for comparison to expected theoret- 

ical R band flux from the RS. Note that 990123 was an exceptional 
burst-a very small fraction of the parameter space produces op- 
tical flashes this bright. 
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ship is due to the estimations made for r ej and more impor- 
tantly 7i r in the analytical expression. 7^ is more difficult 
to estimate in the non spreading case, since the range of 
RS strengths ranges from ultra-relativistic to Newtonian. 
Since we use a mildly relativistic RS estimate for the ana- 
lytical expression, the analytical expression underestimates 
the ultra-relativistic RS flux and overestimates the Newto- 
nian RS flux. Likewise, in the spreading case, we made an 
estimate of 7^ being mildly relativistic; all spreading cases 
are mildly relativistic to Newtonian, so the scatter around 
a linear relationship is much smaller than for non spreading 
ejecta. There is still a small bit of scatter due to the value 
of yi r that we chose (see Table 1), and the analytical value 
is typically within a factor of a few to 10. 

The analytical expression given in Table 1 is a good es- 
timation of the R band flux at deceleration for both spread- 
ing cases. It is less useful for the non spreading cases; in the 
s = 2 non-spreading case, the ordering of break frequencies 
assumed in the analytical expression is often not applica- 
ble, and the linear relationship does not hold over the entire 
parameter space explored-this occurs in the top right hand 
corner of the plot for this case. The cooling frequency is typ- 
ically lower in s = 2 compared to s = 0, and drops below the 
other two break frequencies frequently in this calculation. 

In Figure HJ we have done the same calculation as 
described above for Figure |S| only here we have plot- 
ted the numerical observer frame R band flux against 
only the parameter dependences from the analytical ex- 
pression for R band flux. For example, for s — 0, A = 
R+ 1 (2rj)) , we have plotted numerically calculated RS 
flux /„ against the analytical combination of parameters 
e^r i e^ +1)/4 L 2 -^n ( Pr+1)/4 r p '-Es2. The scatter (width) of 
the spreading regions on this plot are mainly due to the vari- 
ation of p r , while the scatter of the non spreading regions 
is due to p r and the dynamics parameters. The s = 2 non 
spreading case is again not quite linear, for reasons discussed 
above for Figure |5] 

This plot can be used as a tool to constrain the burst 
parameters using observed R-band flux or upper limits at 
or near deceleration using the more accurate results of the 
numerical calculation. Also, this plot can be used to view the 
ranges of expected R-band flux at deceleration if the early 
afterglow is caused by the RS with baryonic or pair enriched 
ejecta. Lines are drawn onto the plot to compare to the early 
afterglow detection of GRB 990123 at the peak (scaled from 
z = 1.61 to z = 1) and upper limits from ROTSE and Swift. 

We find that optical flash emission from GRB 990123 
falls at the very bright end of the distribution of RS R band 
flux at deceleration for all 4 cases shown in Figure |S] There 
are a great many synthetic bursts in our calculation, espe- 
cially in the case of s — 0, spreading ejecta, with intrinsically 
low R band RS flux levels for burst parameters consistent 
with those found from afterglow modeling. This may indeed 
be the reason that the growing number of GRBs with rapid 
follow up do not show a bright optical flash. 

As an example of how to use Figure|S|to constrain burst 
parameters with optical data or upper limits near deceler- 
ation, we look at the case of GRB 050219a. GRB050219a 
was a 23.6 s long SW/S-detected burst of with a fluence 
of (5.2 ± 0.4) x 10~ 6 erg cm" 2 an d was detected in the x- 
ray by XRT jRomano et al. l l2005). An optical upper limit 
from Swift of V = 20.7 was found 96 s after the burst 



JSchadvet al. 1120051) . We find that the upper limit for this 
burst falls near the bottom of the distribution of the non- 
spreading cases (although the distribution can be extended 
to lower values of using even lower limits for the ranges 
on no /A, and EBr) but is in the middle of the distribution 
for the two spreading cases. It seems that RS flux fainter 
than this upper limit is fairly typical in the case of spread- 
ing ejecta. 

Using Figure HJ we can make constraints on the burst 
parameters of 050219a (assuming that the upper limit of 
V — 20.7 applies at deceleration for this burst and there 
was no extinction). In the s = 2 cases, we can make fairly 
severe constraints on ebt and A*. For example, using the 
plot for the s = 2 non spreadin g case, if we assume z = 1 
(although see lBerger et alJl2005l for the redshift distribution 
for Swift GRBs so far) and p r ~ 2.2 we find (see Table 1) 

e^e^^Al^r Elites < 1.2 x 10" 8 . (19) 

Using the observed burst duration and using the assumption 
2 = 1, tdur = 11.8 s, and estimating E52 from the fluence, 
we set E52 ~ 2. We set Ei r — 0.1, a typical value, as ind i- 
cated from afterglow modeling jPanaitescu fc KumaJ2o'02T) . 
Inserting these values and setting To = 100ro,2, we have 

e^L- ' 3 ^- 1 ^ < 5.3 x 10- 7 . (20) 

Since the dependence of this relation on the number of pairs 
and To, 2 is small, we can say that within a factor of a 
few, Al A e° B i<5.3 x 10~ 7 , implying very small values for A* 
and/or ebt- 

If we set ebt ~ 10 -5 , on the lower end of the 
distribution of values fou nd from afterglow modeling 
JPanaitescu fc Kumad 12002^ . A,<10~ 2 . Also, if we choose 
to make Ei r smaller for this burst than the typical value we 
have chosen, say £; r ~ 0.01, then if A* ~ 1, e_g r <10 -7 , still 
a rather low value. 

If we try the s = 2 spreading case, we find similarly 
severe requirements for A* and ebt- The expression for the 
spreading case, however, is much more sensitive to Fo.2, so 
the limits on A* and ebt are not as robust as in the non 
spreading case above. For A, ~ 1 and To, 2 ~ 1, es r ^10 -9 . 
If we require ebt ~ 10 , then A* ~ 0.1. In summary, for 
the s = 2 density profile to apply, with either spreading or 
non spreading ejecta, the upper limit at 96 s requires low 
values of A, or ebt near deceleration. For the s = cases, 
using the typical values above (£i r ~ 0.1), we find the optical 
upper limit at 96 s to be less constraining on no and ebt- 



4 DISCUSSION 

In this paper, we have calculated RS emission from purely 
baryonic and pair-enriched ejecta. We take into account the 
mildly relativistic nature of the RS (we allow for the full 
range of RS speeds, from Newtonian to ultra-relativistic), 
and self consistently calculate synchrotron self absorption 
and inverse Compton cooling in the ejecta. Additionally, we 
have allowed for the ejecta to be cooled by FS synchrotron 
flux incident on the ejecta and for the low-energy RS photons 
to be absorbed in the FS material as they move outward 
from the source. 

We find that the flux in the R band at deceleration de- 
pends very weakly on the ejecta pair content. The RS radio 
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emission is affected by ejecta pairs the most, however the 
shape of the radio fight curve is not affected. The effect of 
pairs on any singie wavefength fight curve (a reduction in 
the flux with increasing number of pairs per proton) can be 
repiicated by varying the shock microphysics parameters in 
the RS; it is a degenerate probiem. It is impossible to deter- 
mine the pair content of the ejecta from a singie wavefength 
light curve. 

It may be possible to determine if iV±>100 by using ob- 
servations of the RS optical and radio light curve peaks. By 
calculating \ ( see Equation ll7H using observations, one may 
be able to determine if 7V±>100; one cannot determine the 
precise number of pairs per proton with x, but whether there 
are a significant number of pairs present in the ejecta. The 
vaiue of x f° r purefy baryonic ejecta is fairiy constant (for a 
given vaiue of p r ) over a wide range of RS strengths; effects 
that may increase the spread in \ (reducing the effectiveness 
of x as a tool for determining ejecta pair content) include 
absorption of radio RS photons in the FS, contribution of 
the FS emission to the RS flux at the light curve peaks, 
and scintillation in the radio light curve. External (FS syn- 
chrotron) flux incident on the ejecta has little effect on the 
value of X- By using radio observations at higher frequen- 
cies, e.g. at 50 GHz or 250 GHZ (mm), effects of scintillation 
and FS contribution to radio emission can be reduced, and 
X, and therefore the pair content of the ejecta, can be more 
accurately determined. 

iLi et al.l (120031) have also looked at the emission from 
pair enriched ejecta; we agree that it is possible for a large 
number of pairs to be present in the ejecta and that the re- 
sulting R band flux is largely insensitive to the number of 
pairs in the ejecta, however we come to different conclu sions 
regarding the resulting ejecta emission. ILi et alJ l)2003f) pre- 
dict a strong flash in the IR band; we find that the peak 
frequency of the RS is largely unaffected by the pair con- 
tent of the ejecta. Indeed, the injection frequency is reduced 
greatly by even a modest number of pairs, however we find 
that in the most frequent arrangement of frequencies, the 
peak of f v is at v a r, which is not affected by pairs much 
at all. Most of the difference between our calculations re- 
sult from the treatment of RS strength; we have taken into 
accou nt the mi l dly re lativistic nature of the reverse shock, 
where ILi et al.l i2003T) have approximated the RS as being 
highly relativistic. Another contributing factor to the differ- 
ence in our results is the more careful and self consistent 
calculation of synchrotron cooling and self absorption fre- 
quencies that we have done here. 

In conclusion, RS emission from pair-enriched ejecta 
looks very similar to that from purely baryonic ejecta. The 
most promising tool to determine if GRB ejecta is pair- 
enriched, by the measurement of the parameter x defined 
in Equation 1171 depends greatly on well sampled observa- 
tions near the peak of the light curves in the optical and 
the radio; if the optical flash is seen quickly after the burst 
with Swift or other ground based telescopes, radio follow 
up in frequency bands between 50 GHz and 250 GHz (mm) 
from 15 minutes to ~ 1 day should provide some information 
about the pair enrichment of the ejecta. 

This research was supported in part by grants from NSF 
(AST-0406878) and NASA-Swift (NNG05G185) (PK) and 
the US-Israel BSF (TP). 
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APPENDIX A: R BAND FLUX AT 
DECELERATION 

Here we write the full analytic expressions for the R band 
flux at deceleration (in cgs units) for each of the four cases 
in Figure 5. 
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